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ABSTRACT 

The magnetic structure of neutron vortices in the superfluid cores of neutron stars 
is determined assuming that the proton condensate forms a type I superconductor. 
It is shown that the entrainment currents induced by the neutron vortex circulation 
cause the proton superconductor to break into successive domains of normal and 
superconducting regions. The Gibbs free-energy is found in the case in which the 
normal domains form cylindrical tubes coaxial with the neutron vortex. The minimum 
of the energy functional corresponds to a tube radius a ~ 0.1 — 0.5 b, where b is the 
outer radius of the neutron vortex. The magnetic field within the tube is of the order 
of 5 x 10 14 G. 
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1 INTRODUCTION 

The knowledge of the properties of neutron vortex lattice in 
the quantum liquid cores of a pulsar (in the simplest case a 
mixture of neutrons, protons, and electrons, called further n- 
p-e phase) is important for the understanding of the coupling 
mechanisms between the superfluid and normal component 
of the star. The problem differs from the case of vortices in 
a rotating uncharged superfluid mainly due to the fact that 
above the nuclear saturation density, p sat ~ 2.6 x 10 14 g 
cm -3 , the 3 P2 superfluid neutron vortex lattice is embedded 
in a 1 So proton superconductor. Under these conditions, the 
relation between the superfluid mass current, p, and super- 
fluid velocities, v, in the hydrodynamic equations assumes a 
tensor character with respect to the isospin indecies (1, 2): 



Pi 
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(1) 



where pik is the superfluid density tensor, which replaces 
the scalar superfluid density in the Landau two-fluid model. 
Equation (^) shows the entrainment effect, i.e. the transport 
of the mass of both condensates by the superfluid motion 
of each condensate, which arises due to the interaction be- 
tween the quasiparticles with different isospin projections 
(Andreev & Bashkin 1975, Vardanian & Sedrakian 1981). 
In the neutron star cores the entrainment effect leads to 
the appearance of electric currents of protons driven by the 
neutron circulation, which become a source of internal ax- 
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ially symmetrical magnetic fields in pulsars (see Sedrakian 
& Shahabasian 1991 and references therein). 

Except for the early stages of pulsar evolution, the tem- 
peratures of the interiors of the stars, T ~ 10 8 K, are much 
lower than the critical temperatures for the superfluid phase 
transition (T c ~ 10 9 — 10 10 K), therefore the variation of the 
proton gap Ai(T = 0) is along the density profile of the star. 
For typical parameters of the outer core of a a neutron star, 
the protons are type II superconductor (Baym, Pethick and 
Pines 1969). With increasing matter density, the attractive 
interaction between protons in the So partial wave chan- 
nel reduces, leading to a reducton of the proton gap. If the 
gap becomes sufficiently small, the coherence length of the 
proton condensate, £, can become larger than the magnetic 
field penetration depth A (more precisely A/£ < l/v2) and 
the proton condensate will undergo a transition from a type 
II to a type I superconducting state. 

The aim of this paper is the determination of the mag- 
netic field distribution within a single neutron vortex line in 
the case when the proton condensate forms a type I super- 
conductor. The present problem bears a certain resemblance 
to the problem of the explanation of the anomalously large 
magnetic flux observed in the thermo-electrical experiments 
on the hollow superconducting cylinders. Our calculation of 
the Gibbs free-energy parallels developments by Aroutyu- 
nian & Zharkov (1981) and Ginzburg & Zharkov (1993) for 
a hollow superconducting sample in the presence of normal 
currents driven around a cylindrical cavity by the tempera- 
ture gradients. 

The paper is organized as follows. In Section 2 we set 
up the basic system of equations of macroscopic superfluid 
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magnetohydrodynamics for the mixture of neutron and pro- 
ton condensates. In Section 3 these equations are applied 
to a single neutron vortex and the Gibbs free-energy is ob- 
tained for this case. Section 4 shows that the minimum of 
the Gibbs potential corresponds to a neutron vortex with a 
coaxial tube of normal protons confined within a finite ra- 
dius a. The dependence of the radius of the tube and its 
homogeneous magnetic field on the density and the values 
of microscopic parameters is determined. Section 5 contains 
several concluding remarks. 



2 THE GIBBS POTENTIAL OF A 

SUPERFLUID NEUTRON-PROTON 
MIXTURE IN A ROTATING SYSTEM 

Consider the system of rotating superfluid neutron-proton 
mixture with a charge neutralizing background of relativis- 
tic electrons and a vanishing number of quasiparticle excita- 
tions (zero temperature limit) . The kinetic part of the energy 
of neutron superfluid, which is by far the most dominant 
part of the energy of the system, is minimized by a lattice of 
neutron vortices. The superfluid executes a course-grained 
rigid body rotation with an angular velocity fi, supported 
by a vortex lattice of density n — 20, /k per unit area, where 
k is the quantum of circulation. The entrainment effect does 
not changes this result to any considerable extent: the cor- 
rection to the mass current of the neutrons is of the order of 
the ratio of proton to neutron density, which is typically of 
the order of 0.01 — 0.05. At the same time, the resulting mag- 
netic energy density is lower by orders of magnitude than 
the kinetic energy density. Thus, the minimization of the to- 
tal energy of the system can be carried out in two steps. In 
the first step, one minimizes the kinetic energy of rotating 
neutron superfluid. Once the neutron superflow pattern is 
fixed by minimization of the rotational energy, in the sec- 
ond step the thermodynamic corrections to the energy of 
the system are found by the minimization of the Gibbs free- 
energy associated with the unentrained part of the proton 
condensate (the motion of the entrained part follows that 
of neutron superfluid and is determined at the first step of 
the minimization procedure). We shall concentrate on the 
second step of the minimization problem under assumption 
that the proton liquid forms a type I superconductor. A re- 
cent account of the case of type II superconductivity is given 
in Sedrakian & Sedrakian (1995). 

Taking account of the foregoing discussions, the ther- 
modynamical Gibbs potential for a mixture of superfluid 
neutrons and superconducting protons can be written as 



M 



T- - I A- j 12 dV -ft x M, 



(2) 



where j 12 is the entrainment current defined below (equa- 
tion 8), A is the vector potential of the magnetic field, ft 
is the angular velocity of the normal component, and the 
free-energy and the angular momentum of the system are 

T = + ^ y {puv't + 2p 12 V! ■ t> 2 + P22V2) dV 

+ i- f B 2 dV-\ f p(n.r) 2 + nxM, (3) 



r x (p 1 +p 2 ) dV, 



(4) 



where p is the total density of baryonic component. The first 
term in equation (Eh, T c , is the superfluid condensation en- 
ergy. The remaining terms are the sum of the kinetic and the 
magnetic energy of the system. Here the gradient invariant 
superfluid velocities are defined relative to their rigid body 
rotation value fi x r 
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(6) 



where mi, 2 denotes the bare mass, xi, 2 the phase of the su- 
perfluid order parameter in the rotating frame; the isospin 
indecies 1 and 2 refer to protons and neutrons respectively. 
As it can be seen from the relation (^), the diagonal elements 
of the superfluid density matrix, pn and P22 are the densities 
of the unentrained parts of the proton and neutron conden- 
sates, respectively, while the off-diagonal elements P12 = P21 
are the densities of entrained parts; these are related to 
the microscopic variables according to relations (Andreev 
& Bashkin 1975; Vardanian & Sedrakian 1981) 

p^ = TTk pl > pll = TTk Pl ' (7) 

where in the mean-field approximation the entrainment co- 
efficient is k = (m* — mi) /ml, and ml and mi are the 
effective and the bare mass of protons, pi is the total den- 
sity of the proton condensate, and the following sum rules 
hold 



p2 — p22 + Pl2\ pi = Pll + Pl2] P = Pl+P2, 



(8) 



where p2 is the total density of the neutron condensate. 

The magnetic induction B — V x A is determined from 
the Maxwell equation, 



4ty , 



V x B = —tin +3 



(9) 



where the proton supercurrent is split into unentrained and 
entrained parts, respectively, 

e 
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(10) 



On the scale of a single neutron vortex, j 12 is the entrain- 
ment current which is generated by the neutron superfluid 
circulation in the neutron vortex, while j u is the Meissner 
current appearing due to the presence of unentrained pro- 
ton condensate (pn 7^ 0). Applying the curl operation to 
equation (^) one finds the London equation for the case of 
type I proton superconductor 

B + A _2 V x V x B = k(j) 5 i2) (r - r ) , (11) 

where 0o = nhc/e is the flux quantum and ro is the radius 
vector defining the position of neutron vortex in the lattice 
plane and A is the magnetic field penetration depth, 

2 2 

^2 = m x c 



Ane 2 pn 

Further, a change of the origin of the free-energy 

T' = T - T C 2 - ~ J P22 v 2 2 dV+^J p (Or) 2 dV, (12) 
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can be made, since the neutron superfluid velocity, Vz, is 
essentially the familiar solution of equation (^) for a single 
neutron vortex, while Tci - the condensation energy of neu- 
trons - is an additive constant. Finally, eliminating the mag- 
netic induction in favor of the vector potential from equation 
(^) and denoting the condensation energy of protons by T c i, 
we find 



ttA 2 ( 



8ttA 2 V2tt 



Vxi 
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-.4] Jl2 + i(VxA) 2 



dV. (13) 



Equations (g), (g), (jll|) and ( |l3| ) form the basic system of 
equations of the problem. 



3 THE DISTRIBUTION OF EQUILIBRIUM 
MAGNETIC FIELD WITHIN A NEUTRON 
VORTEX 

Type I superconductors in the presence of imposed currents 
(or equivalently magnetic field) break into successive do- 
mains of superconducting and normal regions; the geometry 
of the domains is determined by the symmetry of the prob- 
lem. Assume that the proton superconductor goes over to 
the normal state within a cylinder of a radius a, concentric 
with the neutron vortex, forming a tube of a normal mat- 
ter. (The cylindrical geometry of the tube is dictated, in our 
case, by the symmetry of the neutron vortex). Let the vol- 
ume inclosed within the radius a be V< , while that between 
the concentric cylinders a and b, where b is the outer radius 
of the neutron vortex, be V> . 

Because the proton liquid is normal within the vol- 
ume V<, the superconducting currents vanish in that region, 
j 1L = j 12 = 0. The homogeneous magnetic field B = H cm 
(we neglect for simplicity the paramagnetic effects for pro- 
tons in the normal state) is determined, in this region, by the 
continuity condition B z — H cm at the boundary r = a. The 
integration of equation (^|) over the volume V< is straight- 
forward 



r>< n cm tt 

5 = ^r v< - 



(14) 



Decomposing the last term in equation ( |13[ ) in a sum of 
volume and surface integrals and using equation (^|), the 
Gibbs free-energy in the volume V> becomes 



g > = 
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(15) 



where ai is the condensation energy density of supercon- 
ducting protons. 

Let us turn to the determination of non- vanishing com- 
ponents of vectors A, Vxi and j 12 . The cylindrical sym- 
metry of the problem implies that all unknown quantities 
are functions of the distance r from the axis of the neutron 



vortex and the azimuthal angle ip. From equations. (jE|) and 
( p"o| ) the azimuthal components of vectors Vxi and j 12 are 



... y.y , . lib 

K3i2)v = — , (Vxi)i= = — , 



(16) 



where Q = (eh/2mim,2)pi2 and m is an integer number; the 
second relation follows form the quantization of the proton 
supercurrent circulation around the cylindrical tube. The 
non-vanishing azimuthal component of the vector potential 
A v is expressed in terms of B z component of the magnetic 
induction using equation (J^): 

A hcm + ^Ql +x2 dB 1 
e r c r dr 

The boundary conditions for the magnetic inductions are 
B z (a) — H crn and B z (b) = 0. From the first conditon it 
follows, particularly, that the circulation of the vector A on 
the radius a is 



A ■ dl = na H c 



(18) 



and therefore A v (a) = (a/2) H cm , or combining with equa- 
tion @ 



_ m0 o 87rA 2 2A 2 dB z 
7ra z ca z a dr 



(19) 



The determination of A v function is accomplished by finding 
the function B z (r) from equation ([n]); the latter is a source- 
free in the region r > a. The solution with the boundary 
conditions above is 



B.(r) = 
where 
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and K's and 7's are the modified Bessel functions. Using 
this result, equation (^9|) takes the form 



(20) 

(21) 
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where 
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(24) 



Here H e (r) is the magnetic field intensity created by the en- 
trainment currents. The integrations in equation (113) using 
equations ( |lrj| ) and ( p^ ) are straightforward, and the result 
for total Gibbs potential Q = Q K + Q > is 



-aiV> + 
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where q = p — p and 
JV{a) 2A 2 



p = 



V 



a 2 £' 
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(26) 



where po is the limit of p for A — > oo. Note that equation ( |25[ ) 
has the proper Q — > behaviour when A — > oo for a finite 
(6 — a)/b 1, i.e. of the Gibbs potential is measured from 
its value in the normal state (cf. Aroutyunian and Zharkov 
1981). 



4 THE SIZE OF THE NORMAL PROTON 
TUBE WITHIN THE NEUTRON VORTEX 



As may be seen from the expression (|25[), which contains two 
unknown parameters m and a, the increase in the energy of 
the system associated with the magnetic field is minimal if 
one chooses 



m = p 



ira 2 H e (a) 



(27) 



Note that for the present problem (in contrast to the case 
considered in Aroutyunian and Zharkov 1981) there are no 
energetic barriers for the fulfillment of this condition, since 
for the superconducting protons in the n-p-e phase of neu- 
tron stars hysteresis effects are absent. Imposing the condi- 
tion (G7n we find from equations (123) and (Eq), respectively, 



H e (a) = H c 



-ot\ V> + qua 



H 2 

2 -"cm. 

87T ' 



(28) 



Note that the first condition coincides with the requirement 
j 12 = 0. Using the definition H 2 m /8iv = ct\ one finds 



H cm = H t (a) = 2V2^m. 



(29) 



The relation determining the radius a follows by combining 
this expression with equation (pi|): 



b exp 



2v/2^qT 



Ho 



where 



Hn 



2ttA 2 ' 



Oil 



t2/3 



T 4 /3 H 2 



(30) 



(31) 



Here ^(ef) is the density of states of protons on the Fermi- 
surface. 

Table 1 gives the values of microscopic parameters of 
the proton superconductor for several values of total bary- 
onic density. The dependence of the proton gap and the 
effective masses on the proton Fermi wave- vector, Ufi, are 
from Baldo et al (1992). The relation between fc^i and the 
total baryonic density corresponds to the equation of state 
of Wiringa, Fiks & Fabrocini (1988). The last column is the 
radius of the normal region following from equation (^o|) for 
the neutron intervortex distance b = 4 x 10~ 3 cm, corre- 
sponding to the equilibrium rotation frequency of the Vela 
pulsar Q, = 70.6 s _1 . If the proton gap eventually closes with 
increasing density, the radius of the normal region becomes 
b, i.e. the protons go over in the normal state in the whole 
volume of the vortex. 



5 CONCLUDING REMARKS 

The present paper analyses the implications of a transition 
of a proton superconductor from type II to type I super- 
conducting state for the magnetic structure of neutron vor- 
tex lattice in the core of the neutron star, taking into ac- 
count the entrainment effect. The main result is that, owing 
to the entrainment effect, the type I proton superconduc- 
tor breaks into successive normal superconducting regions. 
We analyzed the Gibbs potential in the case in which nor- 
mal regions are of cylindrical shape, which is most naturally 
implied by the symmetry of the problem. We do not ex- 
clude the possibility, however, that the system may support 
a more complicated pattern of normal-superconducting do- 
mains, and those should be examined in future. 

The transition from the type I to the type II supercon- 
ducting state occurs at densities ~ 3 p sa t- These densities 
are indeed achieved in neutron star models with canonical 
masses ~ 1.4 Mq which are based on moderately-soft 
equations of states (e. g. Wiringa et al. 1988). Though the 
density range occupied by type I proton superconductor is 
not as broad as for a type II superconductor, the density 
profile of the star at high densities is commonly quite flat, 
implying that the type I proton superconductor can occupy 
a considerable fraction of the volume of a superconduct- 
ing core. The magnetic field moment, associated with the 
neutron vortex lattice, would be aligned with the rotation 
vector unless the star is a non-axisymmetrical body at the 
instance of the nucleation of superconducting phase. Since 
this would require rotational velocities close to the mass- 
shedding limit, the nucleation of non-axisymmetrical field 
via this mechanism dose not appear to be a relistic possibil- 
ity. The mean magnetic induction of a single neutron vortex 
is (B) ~ H e {a) (a/b) 2 ~ 10 13 - 10 14 G for densities listed in 
the table; it drops, however, to zero in the limit a — > b when 
the entrainment currents vanish. 

Present results on the ground state configuration of the 
type I proton superconductor can be used to address the 
problems of dissipative motion of neutron vortex lattice, re- 
sulting in a effective coupling of the neutron superfiuid to 
the normal electron liquid, heat generation via irreversible 
processes and other related problems. 

Finally, we note that a situation similar to that dis- 
cussed in the present paper can be obtained in the metallic 
type I superconductors of cylindrical shape when a tempera- 
ture gradient is applied between the axis and the boundary 
of the sample. In principle, the effect of transition of the 
central part of the cylinder to the normal state should be 
observable in this type of experiments. 
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Table 1. Parameter values and the size of the normal region in type I proton supercon- 
ductor 



p 


k F1 


va\ /mi 


Ai 




A 


Ho 


a/b 


(xlO 14 g cm- 3 ) 


(fm-l) 




(MeV) 


(fm) 


(fm) 


(xlO 14 G) 




7.91 


0.85 


0.69 


0.3 


54.2 


41.58 


5.9 


0.14 


8.30 


0.88 


0.68 


0.2 


97.0 


39.20 


6.8 


0.32 


8.56 


0.90 


0.68 


0.1 


197.1 


37.87 


7.4 


0.59 



REFERENCES 

Andreev A. F., Bashkin E. P., 1976, Sov. Phys. JETP, 42, 164 
Aroutyunian R. M., Zharkov G. F., 1982, Sov. Phys. JETP, 83, 
1115 

Baldo M., Cugnon J., Lejeune A., Lombardo U., 1992, Nucl. Phys. 
A536, 349 

Baym G., Pethick C, Pines D., 1969, Nat, 224, 673 
Ginzburg V. L., Zharkov G. F., 1993, Journ. Low Temp. Phys., 
92, 25 

Sedrakian A. D., Sedrakian D. M., 1995, Ap. J., 447, 305 
Sedrakian D. M., Shahabasian K. M., 1991, Sov. Phys. Uspechi, 
34, 555 

Vardanian G. A., Sedrakian D. M., 1981, Sov. Phys. JETP, 54, 
919 

Wiringa R. B., Fiks V., Fabrocini A., 1988, Phys. Rev. C38, 1010 



© 1997 RAS, MNRAS 290, 203-207 



